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FOREWORD 


This  paper  examines  a two-sided  strategic  nuclear  exchange 
model  involving  ICBMs,  SLBMs,  bombers  and  societal  value  on 
each  side.  The  first  striker  allocates  his  weapons  against  the 
second  striker's  weapons  and  value  in  such  a way  that,  after 
massive  retaliation  by  the  second  striker,  the  difference  between 
the  second  striker's  value  destroyed  and  that  of  the  first  striker 
is  maximized.  The  paper  presents  computational  results  including 
the  sensitivities  of  the  outcomes  to  various  force  characteristics. 

The  paper  is  a contribution  to  the  literature  in  that  it 
incorporates  realistic  weapon-on-weapon  attrition  functions  into 
a formulation  of  the  two-strike  strategic  exchange  problem  for 
which  globally  optimal  allocations  can  be  guaranteed. 
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Chapter  I 
INTRODUCTION 


Strategic  weapons  planning  models  based  on  two-strike  nu- 
clear exchanges  have  been  an  integral  part  of  strategic  analysis 
since  the  early  1960s.  One  line  of  research  in  this  area  has 
been  the  development  of  models  for  which  one  could  find  "optimal" 
allocations  of  weapons  to  targets.  These  models  frequently  took 
the  form  of  max-min  problems  wherein  the  first  striker  has  a 
choice  of  allocating  his  weapons  against  his  opponent's  value 
targets  (economic  and  industrial  targets),  thereby  obtaining 
immediate  benefit,  or  of  allocating  those  weapons  against  his 
opponent's  weapons,  in  the  hope  of  reducing  the  impact  of  the 
expected  retaliation.  This  question  of  balancing  counterfcrce 
and  countervalue  targeting  was  perhaps  academic  when  strategic 
weapons  were  relatively  invulnerable.  However,  with  improved 
ICBM  accuracy,  this  issue  becomes  increasingly  significant. 

Among  the  approaches  used  in  the  past  to  answer  this 
problem  were  the  Strategic  Weapons  Exchange  Models,  reported 
in  [1]  and  the  Arsenal  Exchange  Model  [2].  Neither  of  these 
models  couM  guarantee  the  globally  optimal  solution  of  the 
max-min  problem  they  considered.  Both  incorporated  Lagrangian 
solution  procedures. 

A branch-and-bound  procedure  was  described  by  Bracken, 

Falk,  and  Miercort  in  [3]  in  which  they  constructed  a piecewise- 
linear  approximation  to  their  original  two-strike  nuclear  ex- 
change model  and  found  the  global  solution  to  the  piecewise- 
linear  model  using  an  algorithm  of  Falk  [5].  The  Bracken,  Falk, 
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Miercort  paper  established  that,  in  principle , this  approach 
could  be  practically  applied  in  two-strike  modelling. 

This  paper  extends  the  Bracken,  Falk,  Miercort  approach 
to  a model  of  more  interesting  detail,  explicitly  modelling 
ICBMs,  submarines,  and  bombers.  Realistic  attrition  equations 
are  developed,  and  the  results  of  numerous  computer  runs  of  the 
model  are  presented  to  demonstrate  some  of  the  remarkable  re- 
sponses of  the  model  to  parameter  excursions. 


I 


I 


I 


Chapter  II 

THEORETICAL  DEVELOPMENT 

Let  us  call  the  first  striker  Red  and  his  opponent  Blue. 

Our  scenario  goes  as  follows:  Red  strikes  first  with  all  his 
weapons  on  Blue's  weapons  and  value,  where  the  term  "value" 
represents  Blue's  capability  to  sustain  its  society,  and  may 
include  population,  industry,  and  so  forth.  Blue  then  retaliates 
with  all  his  remaining  weapons  against  Red  value.  This  inter- 
action is  depicted  in  Figure  1. 

First  Striker  Retaliator 

(Red)  (Blue) 


Figure  1.  STRATEGIC  EXCHANGE 

I 

Red  has  two  concerns:  first,  to  destroy  Blue's  value  and 
second,  to  preserve  his  own  value  through  destruction  of  Blue's 
retaliatory  weapons.  To  investigate  this  further,  we  develop 
• some  notation. 

We  will  denote  Red's  allocation  by  x = (x.^:  i*l,...,I; 

J»1,...,J),  where  x^  is  the  number  of  Red's  class-i  weapons 
directed  against  Blue  class-j  weapons  for  j=l,...,J-l,  and  Xy 
is  the  number  of  Red's  class-i  weapons  directed  against  Blue 
value.  Let  Bj(x)  represent  the  number  of  Blue's  class-j  weapons 
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that  survive  an  attack  x,  for  j=l,...,J-l,  and  let  = (y-^,  • • • 
be  some  Blue  retaliation.  If  Dg(x)  represents  the  damage  to 
Blue  value  resulting  from  attack  x and  DR(^)  is  the  damage  to 
Red  value  resulting  from  Blue  retaliation  then  we  use  as  our 
objective  function  the  difference 

Db(x)  - D R(jr). 


Blue  desires  to  minimize  this  difference  for  any  given  x>  while 
Red  wishes  to  choose  hi s attack  x so  that  this  expression  is 
maximized,  taking  Blue's  retaliation  into  account.  If  m,  , 
i=l,...,I  represents  Red's  initial  inventory  of  class-i  weapons, 
then  this  two-strike  war  can  be  modelled  by  the  following  mathe- 
matical programming  problem 

max  min  Dg(x)  - DR(^) 

V 'r 

subject  to 


J 


l 

j=l 


< 


yj  < BjU) 


xij'yJ  4 0 


1 1 j • • • y I 

j=l,. >.,J-1 

i*l, j =1 > • • • > J • 


This  problem  can  be  greatly  simplified  if  we  make  the 
rather  weak  restriction  that  DR(y_)  be  monotone  nondecreasing 
in  each  component  of  y;  i.e.,  if  Blue  retaliates  with  any  more 
of  any  weapon  type,  then  the  damage  to  Red's  value  does  not 
decrease.  With  this  restriction  the  problem  becomes 

max  Dg(x)  " Dpt(M(x)) 


subject  to 


l| 


J 

l X.  . < m 1=1, ... ,1 

j=l  J 

x.  ■ — 0 1=1, ...  j I;  j = 1 , • • • > J 

B(x)  = {B^x),  B2(x)  , . . . ,BJ_1(x) } . 

This  reformulation  is  crucial  to  our  approach.  The  next  two 
sections  develop  the  functions  (x ) , . . . ,BJ_1  (x ) , DR(^),  anti 
Db(x). 
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Chapter  III 


WEAPON-ON-WEAPON  INTERACTIONS 


This  chapter  develops  formulae  for  attrition  of  Blue  weapon 
resources  resulting  from  a Red  attack. 

For  the  purposes  of  this  paper,  we  will  allow  Red  and  Blue 
four  weapon  classes  each.1  The  weapon  classes,  with  their  dis- 
tinguishing characteristics,  are  as  follows: 


Class  1 — ICBMs.  ICBMs  are  normally  housed  in  silos,  which 
are  attacked  individually.  The  locations  of  the  silos  are 
assumed  to  be  known,  and  they  may  or  may  not  be  defended. 

Class  2--SLBMs-AT-SEA . SLBMs-at-sea  comprise  those  missiles 
on  submarines  that  are  patrolling  on  station.  Our  model 
permits  these  submarines  to  be  attacked  by  strategic  weapons. 
Submarines  are  attacked  individually,  the  attacker  is  un- 
certain of  their  location,  and  they  are  not  defended. 

Class  3--B0MBERS.  Bombers  are  attacked  in  groups,  at  air- 
fields of  known  location,  where  they  may  or  may  not  be 
defended.  (A  "group"  of  weapons  is  a configuration  such 
that  an  attack  against  any  one  of  them  is  an  attack  against 
all  of  them.)  Unlike  missiles,  bombers  may  be  scrambled 
on  warning  and  called  back,  if  necessary.  The  model  assumes 
five  or  fewer  separate  Blue  airfields. 

Class  4--SLBMs- IN-PENS.  This  class  consists  of  SLBMs  in 
submarines  that  are  available  for  strategic  retaliation  but 
are  in  pens  of  known  location.  Thus,  they  must  put  to  sea 
before  they  can  retaliate  effectively.  All  submarines  at 
any  given  pen  are  attacked  as  a group  and  may  or  may  not 
be  defended.  We  allow  for  three  or  fewer  Blue  pens. 


1Problems  involving  more  weapon  classes  can  be  and  have  been  solved.  For 
instance,  the  IDASNEM  model  [4]  developed  at  the  Institute  for  Defense 
Analyses,  using  the  approach  outlined  in  this  paper,  allows  nine  weapon 
classes  per  side. 


If  Red  attacks  Blue,  and  there  are  four  weapon  classes  on 
each  side,  then  there  Is  a maximum  of  sixteen  weapon-on-weapon 
Interactions.  However,  because  one  would  pay  a m pro- 

longed run  times  If  one  were  to  allow  all  Interactions,  it  is 
desirable  to  disallow,  a priori,  certain  unlikely  weapon-on- 
weapon  interactions.  For  instance,  it  is  unlikely  that  Red 
bombers  would  be  used  against  Blue  ICBMs,  since  there  would  be 
sufficient  warning  before  the  bombers'  arrival  to  launch  all 
the  ICBMs,  thereby  removing  them  as  a target.  All  weapons  may, 
of  course,  attack  value  targets. 

We  will  assume  that  Red  ICBMs  and  SLBMs-at-sea  may  attack 
any  Blue  weapon  class.  Red  bombers,  however,  will  only  be 
effective  against  Blue  SLBMs-in-pens , provided  we  assume  that 
the  latter  require  sufficient  time  to  put  to  sea  and  so  are 
vulnerable  to  bombers.  Red  SLBMs-in-pens  may  only  attack  value 
targets.  These  relationships  are  summarized  in  Figure  2. 


x y 


Figure  2.  WEAPON-ON-WEAPON  INTERACTIONS 
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The  following  set  of  equations  gives  the  Blue  surviving 
force  after  a Red  attack.  The  attacking  units  mentioned  below 
are  assumed  to  be  warheads — when  the  attacking  weapon  is  a 
missile--or  bombs  (or  SRAMS,  for  appropriate  parameter  values) 
when  the  attacking  weapon  is  a bomber.  Note  that  we  are  re- 
placing what  are  essentially  random  variables  with  their  expecta- 
tions . 


A.  SURVIVING  BLUE  ICBM  WARHEADS 

n^:  number  of  silos, 

p^:  probability  of  single  shot  kill  of  a silo  due  to  a 
single  attacking  unit  of  class  r, 

p : individual  reliability  of  an  attacking  unit  of 
class  r, 

q : overall  reliability  of  weapon  class  r;  i.e.,  there 
is  a probability  1-q  that  no  weapons  of  class  r 
will  function,  r 

a1,:  probability  that  an  attacking  unit  of  class  r will 

r penetrate  silo  defenses, 

T 

w : number  of  warheads  per  missile, 

x , : number  of  attacking  units  of  class  r assigned  to 

r category  "ICBMs." 

We  require  that  attacking  units  of  any  single  type  be  uni- 
formly distributed  among  the  targets  of  a single  category. 
Therefore,  the  number  of  attacking  units  assigned  to  each  silo 
is  x^/n^.  Then,  the  probability  that  an  individual  silo  sur- 
vives an  attack  by  weapon  class  r,  provided  weapon  class  r does 
not  fail,  is 


/ 

V 


i-p. 


pi\*n/nI 
r ) 


Taking  into  account  the  probability  that  weapon  class  r might 
fail,  the  probability  that  a silo  survives  an  attack  by  weapon 
class  r is 
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(i-qr> 


+ q- 


•K  •; »: ) 


x , /n' 
rl 


Therefore,  recalling  that  only  Red  weapon  classes  1 and  2 will 
attack  Blue  ICBMs,  the  probability  of  a silo  surviving  the  full 
attack  is 


jjd-qp  +ar(i-»r*JpJ)Xrl/n  | 

and  so  the  approximate  expected  number  of  surviving  Blue  ICBM 
warheads  is 

X * vi)WnI\- 


B.  SURVIVING  BLUE  SLBM-AT-SEA  WARHEADS 

3 

n : number  of  submarines-at-sea 

xS:  probability  that  Red  will  correctly  acquire  the 

location  of  an  individual  submarine, 

3 

p : probability  of  single  shot  kill  of  a submarine 
due  to  a single  attacking  unit  of  class  r, 

p : individual  reliability  of  an  attacking  unit  of 
r class  r, 

q^:  overall  reliability  of  class  r, 

3 

w : number  of  warheads  per  missile, 

3 

Y : number  of  missiles  per  submarine, 

x ,5:  number  of  attacking  units  of  class  r assigned  to 
r''  the  category  "SLBMs-at-sea.  " 

As  with  ICBMs,  we  assume  that  the  number  of  class-r  attack- 

3 

ing  units  assigned  to  each  submarine  is  x^/n  . Considering 
only  class-r  attacking  units,  the  probability  that  a single 
submarine  will  survive  is 
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K p?)Xr2/n  . 

if  weapon  class  r does  not  fail.  Taking  this  possibility  into 
account,  the  probability  of  submarine  survival  is 


d-qr)  + Qr  (l-Pr  Pr 


S\  xr2/n 


Now,  the  probability  that  the  submarine  will  survive  an  attack 
by  Red  weapon  classes  1 and  2 is 


n (i-qj  + qji-p„  p 


s\xr2/n‘ 


r r 


All  the  above  assumes  that  the  location  of  the  submarine  is 
known.  Considering  the  possibility  that  the  submarine's  location 
will  not  be  acquired,  the  probability  of  submarine  survival  is 


(1-t  ) + t n (1-q  ) + q„ ( 1_P„  P 


S \ xr2/n 


r r 


Hence  the  approximate  number  of  surviving  Blue  SLBM-at-sea 
warheads  is 


nSwV  (1-tS)W 


2 . 
n 4(i- 

’=i  1 


q ) + q { 1— p p 
Hr'  HrV  ^r  ^r 


C.  SURVIVING  FLUE  BOMBS 

T: 

r r : number  of  bombers  at  airfield  n (n=l,...,5)  (this 
may  be  0 when  there  are  fewer  than  5 airfields), 

n 

n : probability  that  a bomber  will  successfully 
n scramble  from  airfield  n, 

B 

p : probability  that  a bomber  at  airfield  n will  be 
n*  destroyed  when  airfield  n is  attacked  by  a single 
class-r  attacking  unit, 

p^:  individual  reliability  of  a class-r  attacking  unit, 
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q^:  overall  reliability  of  weapon  class  r, 

n 

a : probability  that  a class-r  attacking  unit  will 
5 penetrate  the  defenses  of  airfield  n, 

g 

Yn:  number  of  bombs  per  bomber  at  airfield  n, 

x number  of  class-r  weapons  allocated  to  "bombers" 

r3  (r=l, 2 ) . 

We  will  require  that  attacking  units  be  allocated  propor- 
tionally to  the  number  of  bombers  at  each  airfield,  thus  airfield 
n receives 

mn  x V 1 

n r3  nil  n 

weapons  of  class  r.  The  probability  that  a bomber  at  airfield 
n survives  a class-r  attack,  assuming  no  scramble  and  no  overall 
failure  of  weapon  class  r,  therefore  is 


K c) 


n r3  nil  n 


Including  the  probability  of  successful  scramble,  the  probability 
of  survival  is 

B „ . i B 

B , n B-,  /,  B B \ mn  Xr 3/ nL  mn 
nn  + (1-nn)  \1_pr  °n,r  Pn,rj  n_1 

Recalling  that  the  probability  that  weapon  class  r will  fail  as 
a class  is  1-q^,  the  probability  of  survival  of  a single  bomber 


Hence,  the  probability  that  a single  bomber  will  survive  the 


full  attack  Is 
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2 

n 

r=l 


d-qr) 


+ q> 


(nB  + (i-nB>  (i-p  oB  PB 

\ n n y r n,r  *n,r ) ( j 


The  approximate  number  of  bombs  surviving  at  airbase  n is 

B B 

Y°m^  times  the  above  expression,  so  that  the  total  number  of 
bombs  remaining  is 


2 

n 

r=l 


U-V  + ‘3r{'V  + (l-$  (l-Pr^>r 


D.  SURVIVING  BLUE  SLBM- IN  -PENS  WARHEADS 

D 

m : number  of  submarines  in  pen  n (n=l,2,3)  (may  be 
0 when  there  are  fewer  than  3 pens), 

P 

p : probability  of  single  shot  kill  of  a submarine 
n*  in  pen  n due  to  a class-r  attacking  unit, 

p : individual  reliability  of  a class-r  attacking 
unit , 

q^:  overall  reliability  of  class  r, 

P 

a : probability  that  a class-r  attacking  unit  will 
n,r  penetrate  the  defenses  of  pen  n, 

P 

w^:  number  of  warheads  per  missile  in  pen  n, 

P 

Yn:  number  of  missiles  per  submarine  in  pen  n, 

x i,:  number  of  class-r  attacking  units  assigned  to 
"SLBMs-in-pens"  (r=l,2,3)- 

As  with  bombers,  pens  will  be  attacked  in  proportion  to 
the  number  of  submarines  in  each.  Thus,  the  number  of  class-r 
attacking  units  assigned  to  pen  n is 

P 3 p 

”n  W J mn  • 
n=l 

If  class  r does  not  fail,  the  probability  that  a submarine  in 
pen  n will  survive  all  class-r  attacking  units  is 
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C 


(l-p  aP  pP  ) 
^ Hr  n,r  ^n,r  ) 


m 


P , r P 
n mn 


Including  the  overall  reliability,  the  probability  that  a sub- 
marine will  survive  is 


( 


p . i p 
m x_u/  l - 


/ P P \mn  Xr4'  ^ mn 

1-q  ) + q (l-p  a pr  1 n 4 n=l  n 
4r ' Mr  Hr  ntr  *n,r  ) 


Thus,  the  probability  of  the  submarine  surviving  the  full  attack 
is 


3 

l 

r=l 


l { (l-qr)  + 3r  ^1_Pr  an,r  Pn,r^ 


mn  WJ=1  % 


so  that  the  approximate  number  of  SLBM-in-pens  warheads  that 
survive  is 


jx  “nYnmn  j U'V  + qr  ( °n  ,r  Pn,r ) " ri|  n'1  "J- 
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Chapter  IV 

VALUE  DAMAGE  FUNCTIONS 


For  the  purposes  of  this  paper,  let  us  suppose  that  damage 
to  value  on  both  sides  can  be  adequately  described  by  the 
following  two  functions: 


db(x)  = VB 


(1"exp('J1 


Dp(Z)  * V. 


f]  -exp  (-  l y 

V VJ=1 


J yJ 


where  Vn  and  V0  are  the  aggregate  Blue  and  Red  values  and  the 
^ B n R B R 

parameters  ui,  y ^ and  Vj  have  been  determined.  This  deter- 

mination can  be  made  using  curve  fitting  techniques  on  the 
results  of  a one-sided  allocation  model  for  samples  of  possible 
Red  attacks  and  Blue  retaliations.  It  is  not  even  a requirement 
of  this  approach  that  these  functions  be  expressible  in  explicit 
form.  For  instance,  the  IDASNEM  model  uses  the  output  of  a sub- 
routine that  assesses  the  damage  against  up  to  two  thousand  turret 
classes  per  side. 


Chapter  V 
SOLVING  THE  MODEL 


The  foregoing  allows  us  to  write  explicitly  the  model: 


maximize  V. 


B 


(i-exp ("i^i  u?  xi5  1))-  vR(1-exp(-JI1 


„R  y 

Vi  ^ 


)) 


subject  to 
5 


l x..  < m.  i=l, ... 

1=1  1 


II  ^ j / I I \ xil//n 

yx  = w n n jd-q1)  + q^i^  °±  pJ  j 


y0  = n w Y (1-T  ) + T 


. „SJSVS  x TS 


{(^i5  + qi(i-Pi  pi ) l2/  }j 


' l ^ L K>  * ^ {'fi  + ^ K <i  ^ ^ 


3 

yjt a I 
n=l 

P P P 

WVi 

Xy  >° 

i-l». 

IV 

o 

J-l,. 

we 

solve 

■£  W I.  -P 


we  first  Introduce  new  variables  and  perform  some  elementary 


manipulations  to  put  the  problem  in  separable  form — in  which  each 
function  g(£),  where  ± 5 ),  is  of  the  form 

m 

S(£)  - l SiCq)- 
i=l 

In  separable  form,  this  problem  is 

maximize  ]vn (l-exp (-£B) j - Vn  ( 1-exp  (-E;R) )[ 
x z ^ ' RV 


subject  to 


4 „ v 


B 


< l u!  x,  1 


rR  > T uR  z VJ 
? yJ  Zj 


1=1 

4 


■i  Ai5 


R 


5 

l xiJ  s mi 


i-1, . . . ,4 


In  z1  - ln(wInI ) > ^ In  | ( l-q± ) + Q^l-P.^  0i  pJ  ^ 11  j 


“[(*•“•■■■)?]*,!,  “K>- I 


■j  > j,  >: «:  'i 


In  c?  > I In 
n i-1 


a-%) + % K + ^n> K <i  <i I""  } ( 


n»l,...,5 


zn  > l wP  yP  mP  CP 
*♦  - n 'n  n n 


16 


P 3 p 

p 3 ( / P P \m  - xi h/  £ % i 

lnCn-  ta|(1-qi)+qi(1-pi0n,iPn,i)  ' k=1  j tFl--’3 

R B B P 

All  , z y E,  , £ , C“,  > 0.  Note  that  we  have  changed  some 

equalities  to  inequalities.  This  was  done  only  to  improve  the 
efficiency  of  the  computation.  It  is  apparent  by  monotonicity 
that  all  inequalities  will  be  equalities  at  the  optimal  solution. 

The  next  step  is  to  place  upper  and  lower  bounds  on  all  the 
variables.  This  can  easily  be  done  because  the  allocations  on 
each  side  are  bounded  by  the  available  resources.  The  bounds 
appear  in  Table  1.  In  ordc  - to  save  space  in  this  table,  we 
have  occasionally  employed  the  notation  £(c)  and  u(c)  to  denote 
the  lower  and  upper  bounds  on  variable  c,  respectively. 

Because  the  problem  comprises  only  scalar  functions  of 
bounded  variables,  and  because  these  functions  remain  bounded 
over  their  domains,  we  can  approximate  these  functions  by 
piecewise-linear  functions  in  the  following  way. 

Each  function  g(?)  in  the  separated  problem  is  a function 
of  a scalar  variable  for  which  we  have  lower  and  upper  bounds 
£(?)>  u(c).  We  divide  the  interval  [£(c)»  u ( c ) 3 into  subintervals 
determir^d  by  the  "cut  points’  t°,  tl,...,t*  where  t°  = i.(C)» 
t^  = u(c)  and  t°  < t1  < t2  < . . . < t^.  Then  we  compute  g(t^) 
for  i = 0,...,k  and  form  the  piecewise-linear  function 

g(t>  * -rrr^r  [s<tm>  - g(t1)]  + set1) 

for  t E [t1.  ti+1]  i-0 k-1. 

This  process  is  depicted  in  Figure  3* 

The  problem  thus  generated  will  be  called  the  "piecewise- 
linear  two-strike  problem."  Since  all  functions  in  the  separated 
problem  are  continuous,  so  are  all  the  functions  of  the  piecewise- 
linear  problem.  Further,  by  taking  enough  cut  points  for  each 
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Table  1.  LOWER  AND  UPPER  BOUNDS 


variable,  the  piecewise-linear  two-strike  problem  can  be  made  to 
approximate  the  original  two-strike  problem  arbitrarily  closely 
1 according  to  any  of  a number  of  standard  measures  of  closeness. 

The  piecewise-linear  problem  has  the  advantage  that  it  can 
be  solved  using  a practical  branch-and-bound  algorithm  developed 
by  Falk  [5].  This  algorithm  has  been  programmed  for  the  computer 
twice.  The  NUGLOBAL  code  [7]  used  by  Bracken,  Falk,  and  Miercort 
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tk=  u(C) 


Figure  3.  PIECEWISE-LINEAR  APPROXIMATION 

I 

developed  problems  when  applied  to  the  model  we  have  described 
above.  A more  stable  code  was  developed  by  Grotte  [6]  and  was 
the  code  used  in  the  runs  described  below. 
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Chapter  VI 
RESULTS 

The  piecewise-linear  model  was  run  on  a CDC  6^00  computer. 

For  all  runs,  upper  and  lower  bounds  were  computed  as  described 
above  and  eight  cut  points  per  variable  were  chosen  as  follows. 
Noting  that  all  functions  g(c)  in  the  separated  two-strike  model 
had  derivatives  that  were  largest  in  absolute  value  for  small 
C,  and  that  decreased  in  absolute  value  as  c increased,  wc  felt 
that  small  intervals  for  small  c , and  large  intervals  for  large 
C,  should  be  used  for  constructing  the  piecewise-linear  approxi- 
mating function  g(c).  Therefore,  4(c)  and  u(c)  were  computed 
first,  r = u(c)  - (c)/  127  was  then  calculated,  and  the  cut 
points  for  c were  chosen  to  be 

t1  = 4(c)  + (21-l)r  i-0,  ... ,7- 

In  this  way,  each  subinterval  was  twice  the  length  of  the  one 
preceding  it,  ensuring  that  there  would  be  small  sub  intervals 
at  the  low  end  of  the  interval  [4(c),  u(c)]  and  larger  sub- 
intervals at  the  high  end.  In  practice,  this  procedure  produced 
better  results  than  using  many  more  intervals  of  uniform  length. 

Our  investigation  of  the  responses  of  the  model  begins  with 
a base  case  set  of  parameters.  These  are  listed  in  Table  2.  A 
number  of  series  of  runs  were  performed  in  each  of  which  a specific 
subset  of  parameters  was  varied  in  a systematic  fashion.  Certain 
model  outputs  were  then  graphed  to  demonstrate  the  sensitivity 
of  the  model  to  parameter  changes.  These  are  discussed  below. 

Each  run  took  between  100  and  200  seconds  of  CPU  time. 
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Table  2.  BASE  PARAMETERS 


VR 

= 

5000.0 

m1  = 6000.0 

qx  = 0. 

9999 

pl 

= 0.8 

VB 

= 

5000.0 

m2  = 1800.0 

9 2 “ 0 . 

9999 

P2 

= 0.8 

= 1500.0 

^3  = 0 • 

9999 

p3 

= 0.9 

m4  = 1200.0 

q4  ■ o. 

9900 

P/. 

= 0.7 

M? 

= 

0.0005 

yj  = 0.0005 

Pi  = 0. 

7 

-! 

= 0.7 

y2 

= 

0.0003 

y2  = 0.0003 

B _ n 
y 2 = 0. 

5 

p2 

= 0.5 

B 

p3 

= 

0.0003 

y!*  = 0.0006 

U®  . 0. 

5 

R 

p3 

= 0.5 

J 

4 

= 

0.0003 

yj  = 0.0003 

yB  = 0 

5 

y4 

= 0.5 

n1 

= 

1000.0 

pj  « 0.8 

o|  = 0. 

95 

w1 

= 

6.0 

= 0.7 

al2  = 0.95 

nS 

= 

50.3 

Pi  = 0.9 

wS 

= 

10.0 

P?  = 0.9 

s 

Y 

= 

16.0 

s 

T 

= 

0.3 

N8 

= 

5.0 

yB 

Y1 

= 

10.0 

"l 

= 100.0  n i = 

0.5 

0n  ■ 

0.95 

»!*■ 

0.95 

Y* 

= 

10.0 

mB 

= 100.0  n2  = 

0.  5 

8Si  ■ 

0.95 

a22  = 

0.95 

YBj 

= 

10.0 

m3 

= 200.0  nB  = 

0.4 

°31 

0.95 

°33 

0.95 

yb 

t4 

= 

10.0 

m4 

= 200.0  n4  = 

0.4 

°41 

0.95 

aB  = 

a44 

0.95 

yB 

= 

10.0 

m5 

= 50.0  rig  = 

0.6 

a®  = 
51 

0.95 

_B  r 
a55  * 

0.95 

t 


Table  2.  BASE  PARAMETERS  (continued) 


P?1  - 

0. 

97 

p?2 

= 0, 

,97 

II 

r-H 

CO  CVJ 
Q. 

0. 

,97 

P22 

= 0. 

,97 

p31  = 

0, 

,97 

P 32 

= 0, 

,97 

pSi  - 

0. 

97 

p42 

= 0, 

,97 

Pli  - 

0. 

,97 

= 0, 

,97 

NP  = 3.0 


wP  = 10.0 

4 

= 16.0 

mP  _ 
ml  " 

25.0  aP 

= 0.9 

= 0.9 

•I 

= 0.9 

wP  = 10.0 

= 16.0 

P _ 
2 " 

25.0  aP 

= 0.9 

= 0.9 

= 0.9 

wP  = 10.0 

Y3 

= 16.0 

P 

3 

25.0  aP 

= 0.7 

4 

« 0.7 

4 

= 0.7 

- 0.85 

P12  = 

0.85 

0.85 

'21  = 0.85 

P22 

0.85 

p23  = 

0.85 

'31  = °'95 

P32  = 

0.95 

P33 

0.95 

A.  SERIES  1 --MULTI PLES  OF  RED'S  TOTAL  ARSENAL 

The  first  series  of  runs  demonstrates  the  effects  of  changes 

• in  the  number  of  weapons  in  each  Red  weapon  class.  This  was 
achieved  by  multiplying  the  basic  values  of  the  parameters  m^, 
m^,  m^  and  m^  by  0.2,  0.3,  0.4,  0.5,  0.6,  0.7,  0.8,  0.9,  and 
1.0  so  that  the  range  of  20  percent  to  100  percent  of  the  basic 

* Red  arsenal  was  covered.  Inspection  of  the  resulting  damage 
curves  in  Figure  4 shows  that  as  Red's  arsenal  increases,  Red's 
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Figure  4.  VALUE  DAMAGE  AS  A FUNCTION  OF 
MULTIPLIER  OF  RED  ARSENAL 


ability  to  destroy  Blue’s  value  increases  almost  linearly , but 
Red’s  ability  to  protect  itself  from  Blue  retaliation  does  not 
improve  for  force  levels  above  about  50  percent  of  the  basic 
arsenal  size.  Figure  5 indicates  how  Red  allocates  weapons 
between  counterforce  and  countervalue  missions  for  various 
arsenal  sizes.  Note  that  for  small  arsenals.  Red’s  optimal 
strike  involves  a relatively  high  level  of  counterforce  operations. 


Figure  5.  PERCENT  RED  WEAPONS  ASSIGNED  TO  COUNTERFORCE  MISSIONS 
AS  A FUNCTION  OF  MULTIPLIER  OF  RED  ARSENAL 
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B.  SERIES  2--CHANGES  IN  RED  ICBM  AGAINST  BLUE  SILO  KILL 

PROBABILITY 

For  the  second  series  of  runs,  we  changed  only  p^,  which 
represents  the  probability  of  single  shot  kill  of  a Red  ICBM 
warhead  against  a Blue  silo.  Figure  6 exhibits  the  resulting 
value  damage  curves.  It  is  apparent  that  for  values  of  p^ 
below  about  0.5,  changes  in  ICBM  hard-target  kill  probability 
make  little  difference  in  the  outcome.  Observe  that  as  p^ 
increases,  the  damage  differences  increase,  which  is  expected; 
but  total  damage  to  Blue  value  generally  decreases , since  the 
more  effective  Red  ICBMs  are  diverted  from  the  c 'untervalue 
attack  in  a successful  effort  to  reduce  retaliatory  damage  to 
Red.  Note  the  undulatory  behavior  of  the  curves.  We  can  inter- 
pret Figure  6 as  follows:  For  p^  less  than  0.4,  Red’s  attack 
remains  the  same,  resulting  in  level  damage  curves  for  both 

T 

Blue  and  Red.  For  p^  between  0.4  and  0.6,  Red  can  take  advantage 
of  increased  ICBM  effectiveness  and  so  divert  some  of  these 
weapons  from  countervalue  to  counterforce  missions.  With  fewer 
weapons  against  Blue  value,  the  damage  to  Blue  decreases. 

Damage  to  Red  decreases  because  of  increased  counterforce  attack. 
As  pj  increases  between  0.6  and  0.7,  Red  finds  he  can  accomplish 
the  same  counterforce  goal  with  fewer  ICBMs  and  so  he  reapplies 
the  excess  to  countervalue  assignments,  thereby  increasing  the 
damage  to  Blue  while  keeping  the  retaliatory  damage  to  himself 
constant.  For  p^  between  0.7  and  0.75,  Red  can  cause  more 
counterforce  damage  with  fewer  weapons  so  that  damage  to  himself 
actually  begins  to  drop  while  damage  to  Blue  still  increases. 
Between  p*  = 0.75  and  p*  = 0.9,  Red  again  diverts  weapons  from 
value  targets  to  Blue  ICBMs,  so  damage  to  both  sides  decreases. 

T 

For  p^  above  0.9,  the  curves  begin  to  repeat  the  behavior  seen  for 
p.  greater  than  0.6.  This  behavior  results  from  rapidly  changing 
allocations  as  p^  changes,  and  indeed,  one  should  note  that  in 
the  case  of  a "flat"  global  optimum,  or  in  the  case  of  many 
local  optima,  all  with  values  close  together,  the  allocations 
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Figure  6.  VALUE  DAMAGE  AS  A FUNCTION  OF  THE  PROBABILITY  OF  A RED 
ICBM  KILLING  A BLUE  SILO,  p! 
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may  change  almost  discontinuously  as  certain  parameters  are 
varied,  even  though  the  objective  function  changes  smoothly. 

C.  SERIES  3--CHANGES  IN  RED  SYSTEM  RELIABILITY 

This  series  of  runs  was  conducted  by  assigning  the  values 
0.2,  0.3j  0.4,  0.5,  0.6,  0.7,  0.8,  and  0.9*  one  at  a time,  to 
all  parameters  q^,  q^s  and  q^.  These  parameters  represent 
the  system  reliability  of  Red's  four  weapon  classes  against 
Blue  weapons  only.  The  resulting  damage  curves  are  displayed 
in  Figure  7.  Note  the  swift  decline  in  the  damage  suffered  by 
Red  as  compared  with  the  much  slower  decline  in  the  damage  suf- 
fered by  Blue  as  Red  shifts  weapons  of  increased  reliability  to 
counterforce  operations.  Contrast  the  nonmonotonic  behavior  of 
the  Blue  damage  curve  to  the  strictly  decreasing  Red  damage 
curve  and  the  smooth  rise  of  the  damage  difference  curve.  Figure 
8 shows  how  the  surviving  fraction  of  Blue  warheads  changes  for 
each  class.  Observe  that  for  low  Red  system  reliabilities.  Blue 
weapon  class  1 (ICBMs)  survives  better  than  Blue  weapon  class  3 
(Bombers),  while  for  high  Red  system  reliabilities.  Blue  weapon 
class  3 has  the  edge. 

D.  SERIES  4--CHANGES  IN  BLUE  BOMBER  SCRAMBLE  PROBABILITIES 

This  series  tests  the  effect  of  varying  the  Blue  Bomber 
scramble  probabilities.  The  variation  was  effected  by  assign- 
ing the  values  0.5,  0.6,  0.7,  0.8,  and  0.9,  one  at  a time, 
jointly  to  ri°,  n|»  n^,  and  The  damage  curves  appear  in 

Figure  9*  The  result  is  that,  as  ail  the  increase  Red  is 
not  induced  to  changt  his  attack,  but  Blue  can  inflict  greater 
retaliatory  damage  owing  to  his  improved  bomber  survivability. 

E.  SERIES  5 — MULT  I PLES  OF  BLUE'S  TOTAL  ARSENAL 

The  last  series  examines  what  happens  when  Blue's  arsenal 
is  increased  over  the  basic  values,  by  multiplying  the  basic 
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Figure  7.  VALUE  DAMAGE  AS  A FUNCTION  OF  RED  SYSTEM  RELIABILITY 


values  of  all  the  parameters  n , n , , m , rru,  m?,  m: , 

p p ^ j ■ ^ ^ 

m;,,  and  by  1.0,  1.1,  1.2,  1.3*  1.4,  and  1.5.  Figure  10 
displays  the  damage  curves.  While  an  increased  Blue  arsenal 
improves  Blue's  retaliatory  capability,  it  further  serves  a 
secondary  role  of  drawing  Red's  attack  away  from  Blue's  value 
targets. 


• BLUE  WEAPON  - CLASS  1 

* BLUE  WEAPON  - CLASS  2 
■ BLUE  WEAPON  ♦ CLASS  3 

♦ BLUE  WEAPON  - CLASS  4 
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